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ABSTRACT 

We have used a Fourier Transform Spectrometer on the James Clerk Maxwell Tele- 
scope to study the submillimetre continuum emission from dust in three hot molecular 
cores (HMC). The spectral index f3 of the dust emission for these sources has been 
determined solely within the 30 GHz wide 350 GHz (850 iim) passband to an accuracy 
comparable to those determined through multi-wavelength observations. We find an 
average (3 ~ 1.6, in agreement with spectral indices determined from previous submil- 
limetre observations of these sources and with those determined for HMC in general. 
The largest single source of uncertainty in these results is the variability of the at- 
mosphere at 350 GHz, and with better sky subtraction techniques we show that the 
dust spectral index can clearly be determined within one passband to high accuracy 
with a submillimetre FTS. Using an imaging FTS on SCUBA-2, the next generation 
wide-field submillimetre camera currently under development to replace SCUBA at 
the JCMT in 2006, we calculate that at 350 GHz it will be possible to determine 
to ±0.1 for sources as faint as 400mJy/beam and to ±0.3 for sources as faint as 
140 mJy/beam. 
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1 INTRODUCTION 

Stars form in dense cores of dust and gas within molecular 
clouds, but remain cocooned in the material of their na- 
tal envelopes. Though young stars and protostars can be 
extremely luminous, the enveloping dust absorbs nearly all 
of the ultraviolet and optical starlight and reradiates it at 
longer wavelengths. In the submillimetre regime, the dust 
emission is generally optically thin. Submillimetre observa- 
tions thus sample all the dust emission within the cloud 
along the line of sight, and can be used to de termine the 
colum n density and mass of material present dHildebrandl 
Il983fl . Interpretation of the observations, however, depends 
on the assumed dust grain temperature, composition and 
size distribution. The effects of composition and grain size 
on the dust emission are usually combined into a dust opac- 
ity, re„, which can be reasonably fit empirically with a power 
law dependence on frequency in the submillimetre. 

At submillimetre wavelengths, dust is generally as- 
sumed to radiate as a blackbody at some average temper- 
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ature Td modified by the optical depth, r v oc k„. For opti- 
cally thin emission, we then find S v oc Ku B v (Td). Assuming 
a power law dependence on frequency, the dust opacity k v 
can be parameterized for a given opacity /to at frequency uq 
as k v — ko{v I 'vq) 13 , where 3 is the spectral index of the dust 
emission. For optically thin emission, the dust continuum 
emission then follows the power law: 

S„ocko( — fB v {T d ) (1) 

Through mathematical modelling, k„ and 3 have been cal- 
culated for various dust grain compositions and size dis- 
tributions iDraine fc Leelll984HOssenkopf fc Hennindll99l 
iPollack et alJll9 94h which can be compared with astronom- 
ical observations. Often, however, various combinations of 
the dust parameters can reproduce existing observational 
data with similar accuracy, while dust temperatures and 
densities can vary dramatically within star forming envi- 
ronments. Calculations oft en assume 3 ~ 2 a t these wave- 
lengths, based on work bv iHildebrandl <ll983t) . In practice, 
however, theoretical dust grain models find spectral indices 
anywhere between zero and three (see the summary by Gold- 
smith, Bergin & Lis (1997)), depending on the physical dust 
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Table 1. Positions and peak brightnesses for observed objects 



Source 


Q J2000 


<5j2000 


-^450 




Distance 




h m s 


o / // 


Jy/8.5" beam 


Jy/15" beam 


kpc 


G10.47 


18 08 38.2 


-19 51 50 


260 ± 65 


39.0 ± 5.9 


5.8 


G31.41 


18 47 34.5 


-01 12 43 


130 ± 32 


24.2 ± 3.6 


7.9 


G12.21 


18 12 39.7 


-18 24 20 


121 ± 30 


13.0 ± 2.0 


13.5 



The positions, peak brightnesses at 450 um and 850 |om, and distances to each 
of the hot molecular cor es observed in this stu dy, ordered by Tgso- Positions and 
distances are taken from Hatchc ll et alj l2000l) . while peak brightnesses were cal- 
culated from the archived SCUBA data. 



properties of the model, while observations have found a 
similar range of values in different states of the interstellar 
medium. 

In high mass star for ming regions and the cores 
contai ning warm protostars, IWilliams. Fuller fc Sridharanl 
( 2004) find the average frequency dependence of dust 
emission from more than 60 high mass protostellar ob- 
jects to be characterized by j3 = 0.9 ± 0.4. This is 
compara ble to values observed in less evolved, colder 
objects JGoldsmith. Bergin fc Lis! Il997l: IVisser et al.lll998t 
Hoeerheiide fc Sandelll l2000t iBeuther. Schilke fc Wvrowskil 
2004), but is systematically lower than spectral indices de- 
termined through submillimetre observations of hot molec- 
ular cores (HMC) and ultra-compact HII (UCHII) regions, 
which are likely the next phases in the evolution of high 
mass protostars. In these more evolved cores, values have 
been determined to be near toj|Osorio L _Lizano fc D'Alessiol 
119991) or slightly higher than llHunterlll99gft ~ 2. 

The changing values of j3 in these studies are indicative 
of the evolution of the dust itself in the changing physical 
environments associated with the process of star formation. 
The growth of dust grains, the formation of icy mantles, 
and changes in grain composition are likely all significant 
factors in the variation of the frequency dependence of the 
dust emission in these regions. In order to understand the 
properties and evolution of dust in star forming regions, ob- 
servational methods are required that can discern between 
temperature and density effects and actual changes in the 
dust composition and opacity. 

To determine ko, observations must be calibrated 
against known dust column densities, which generally re- 
quires comparison with observations at shorter wavelengths. 
In contrast, the spectral index j3 of the dust opacity can be 
calculated solely in the submillimetre by comparing obser- 
vations at two widely separated submillimetre wavelengths. 
This method provides great leverage on f3, but is complicated 
by calibration issues from observing in the different pass- 
bands, and has the disadvantage that the value determined 
for P is dependent on the assumed dust temperature. The 
Planck function at 450 u.m is far from the Rayleigh- Jeans 
limit for dust temperatures less than ~ 40 K, and even at 
the higher temperatures expected in high mass star forming 
regions, the change in emission between passbands could be 
due to either small changes in temperature or large changes 
in the dust emissivity. The uncertainties in the determina- 
tion of f3 are thus large enough that the subsequently derived 
physical properties of dusty systems, such as the mass es- 
timates of the material in dense cores in molecular clouds, 



may range over factors of a few. For this reason, it is bene- 
ficial to observe in one passband only. 

One instrument that is well-suited to this problem is a 
Fourier Transform Spectrometer (FTS). FTS have several 
advantages and disadvantages compared with other spec- 
trometer designs. The highest spectral resolution of an FTS 
is determined by the maximum optical path difference be- 
tween the two beams of the interferometer, which is set by 
the length of the translation stage. In practice, space limi- 
tations require that most submillimetre FTS provide low to 
intermediate spectral resolution (R~ 200 - R~ 2000 - 3000 
at 345-375 GHz), excellent for observations of sources with 
wide, bright spectral lines and continuum emission. For these 
studies FTS have a distinct advantage over high spectral res- 
olution but narrow bandwidth heterodyne receivers. The in- 
stantaneous bandwidth of an FTS is inherently broad, and in 
the submillimetre is limited only by the atmospheric trans- 
mission windows, enabling a wider spectral coverage within 
these windows than is possible with the heterodyne receivers 
currently in use. In the past decade, two groups have suc- 
cessfully used FTS to study the submillimetre emission from 
several astronomical sources, such as the Ori on molecular 
cloud iSerabvn fc Weissteinl Il995'l . the Sun iNavlor et all 
2000) , and the a tmospheres of planets jNavlor et alJll994l : 
bavis et al.lll997ft . 

A submillimetre FTS operated at low resolution can 
sample the continuum dust emission from a source with a 
large enough bandwidth that the change in dust emissivity 
with frequency, and thus /3, can be determined solely within 
one passband. At very low temperatures, it is still difficult 
to separate the temperature component of the emission from 
the dust power-law emissivity function, but single bandpass 
observations are able to place limits on the combined dust 
temperature and /3. 

In this study, we present low resolution 350 GHz con- 
tinuum observations of four hot molecular cores performed 
with a submillimetre FTS of the Mach-Zehnder design (MZ 
FTS) developed for use at the JCMT by the FTS group 
at the U niversity of L e thbrid ge, AB, Canada (described in 
detail bv lNavlor et~ail J2003t) 1. 

We describe in the following the observational proce- 
dure and analysis of the data. We compare the dust spectral 
index values determined for our sources with those calcu- 
lated from previous submillimetre measurements, and dis- 
cuss the detrimental effects of the variable submillimetre 
atmosphere on our sensitivity. Finally, we comment on the 
continuum measurement capabilities of an imaging FTS cur- 
rently in preparation for use at the JCMT with SCUBA-2. 
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2 OBSERVATIONS 

2.1 Source selection 

The sources for this pilot study were cho s en fro m a sam- 
ple of HMC from work by lHatchell et"aT1 fcOOOl) for their 
large submilli metre fluxes and lack of extended s tructu re as 
shown bv lHatchell et alJ fcOOOl) and lWalsh et alJ J2003I) . The 
sources, G10.47, G12.21 and G31.41, have been previously 
studied in the submillimetre regime with the Submillime- 
tre Common User Bo l ometer Array (SCU BA) at the JCMT 
lHatchell et all l2000l: IWalsh et all 120031) and in the radio 
continuum (|Wood fc Churchwell l|198flf l. among others). All 
three have also been detected in numerous molecular lines 
and maser emi s sion ( see the detailed source descriptions in 
lHatchell et, alJ feflOfli). Table [J giv es source positions and 
distances from lHatchell et alJ lioOCl) . and 450 urn (660 GHz) 
and 850 \im (350 GHz) peak brightnesses calculated using 
the archived SCUBA data 1 . Each source is associated with 
one or more UC H II regions near the peak of the submil- 
limetre emission. The hot core in G10.47 contains three 
UCHII regions, G10.47-0.03A, B and C, while the other 
sources are coincident with one UCHII region each, G12.21- 
0.10 and G31.41+0.31. These HMC are very peaked in the 
850 \im emission, although G31.41 does contain some more 
extended emission, and the hot cores themselves of these 
sources are much smaller than the JCMT beam with diam- 
eters of only 1-2". C ore mass estimates a re on the order of 
a few thousand M lHatchell et alj|2000l) . A fourth HMC, 
G43.89, was observed by us to only S/N~ 1.5, too low to 
analyse with any degree of certainty, and is therefore omit- 
ted from further discussion. Observations of Mars were also 
performed for flux calibration. 

2.2 Observational method 

The observations were performed during 2003 April 21 - 
28, using the MZ FTS mounted on the right Nasmyth plat- 
form of the JCMT on Mauna Kea, Hawaii. The April 2003 
weather was good to excellent with the zenith opacity mea- 
sured at 225 GHz (1.3 mm), T225, ~ 0.07 for most of the 
observations, giving precipitable water vapour (pwv) levels 
averaging ~ 1.1mm. At 350 GHz, the FTS bandwidth is 
~ 30 GHz. The FTS was operated at low spectral resolu- 
tion (1.5 GHz), resulting in ~ 20 resolution elements across 
the band. The full width half maximum (FWHM) of the 
350 GHz beam at the JCMT is ~ 15", which was matched 
by a cold (4 K) aperture located at the field stop of the detec- 
tor system iNavlor et all 1 1999(1 . The atmospheric emission 
and transmission window at 350 GHz is shown in Figure Q 
Spectra were obtained in groups, with five spectra taken 
of a molecular core immediately followed by five spectra 
of the background sky, +1425" in right ascension from the 
source. The offset in right ascension was calculated to com- 
pensate for the change in airmass of the source due to the 
Earth's rotation during the time taken to obtain five spec- 
tra, allowing the sky measurements to be taken at the same 
airmass as the original source observations. This is crucial 

1 Guest User, Canadian Astronomy Data Centre, which is oper- 
ated by the Dominion Astrophysical Observatory for the National 
Research Council of Canada's Herzberg Institute of Astrophysics 
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Figure 1. The atmospheric radiance (solid line) in units of 
Jy/15" beam and transmission (dash-dotted line) calculated for 
precipitable water vapour levels of 1.0mm. 

for accurate background sky subtraction, as the brightness 
and variability of the atmosphere at 350 GHz is easily the 
largest source of uncertainty in our final results. Five scans 
were obtained approximately every 90 seconds. Each tar- 
get source was observed for between three and four hours, 
with approximately half of that time spent on source and 
half spent on background observations. Due to its bright- 
ness, Mars was observed for a half hour. A pointing check 
performed approximately every 45 minutes resulted in small 
(~ 1 — 2") corrections to the pointing. 

2.3 Data reduction 

Preliminary data reduction was perfor med using processin g 
pipeline software written for the FTS iNavlor et all feooa). 
The data were first manually screened for transients, such as 
cosmic rays, which occur approximately once in every ten in- 
terferograms and produce large intensity spikes at one or two 
data points in the interferograms. These spikes are easily de- 
tected and can be removed within the pipeline. Some spikes 
in the interferograms were spread over multiple data points; 
these interferograms were omitted from the final dataset. 
Only a handful of interferograms were rejected for this rea- 
son. Since the FTS optics and detector have negligible dis- 
persion over the range of interest, a linear phase correction 
was applied to each interferogram. The individual interfero- 
grams were then Fourier transformed. The spectra were then 
coadded into the on- and off-source groups of five, and the 
coadded off-source background scans were subtracted from 
the coadded on-source scans. The coadded, sky-subtracted 
spectra of G10.47, G12.21, G31.41 and Mars are shown in 
Figure H 

Each background-subtracted spectrum was then cor- 
rected for atmospheric tra nsmittance w i thin t he band using 
the technique described bv lDavis et all ll997l) . Updated at- 
mospheric opacities were calculated using the University of 
Lethbridge Transmissio n and Radiance Atmospheric Model 
(ULTRAM; IChapmanl feOOlll ). ULTRAM is a radiative 
transfer atmospheric model which calculates the radiance 
and transmission of the atmosphere through a line- by- line, 
layer-by-layer analysis, and was written with the specific 
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Figure 2. The raw, coadded, sky-subtracted spectra of G10.47, G12.21, G31.41 and Mars before calibration or atmospheric 
transmission correction. Error bars are at the 1 a level. Note that the error bars for Mars lie within the line thickness. 



goal of accurately modelling the atmosphere above Mauna 
Kea. The model parameterizes the atmospheric transmit- 
tance in terms of the pwv and airmass, and includes H2O, 
O3 and O2 as atmospheric absorbers. The airmass of each 
scan was recorded with each interferogram. The opacity at 
225 GHz, T225, was recorded by a dipping radiometer oper- 
ated by the Caltech Submillimeter Observatory (CSO). The 
use of this opacity value in the data analysis introduces some 
uncertainty into the final results, as the CSO radiometer 
operates at a fixed azimuth which rarely corresponds with 
the azimuth of the JCMT observations, and only updates 
approximately every 10-15 minutes, while the pwv levels in 
the atmosphere may vary significantly on shorter timescales. 
This can be seen in Figure |3J where T225 from the CSO is 
plotted against T225 determined from measurements of the 
pwv recorded by the JCMT's water vapour monitor (which 
is pointed slightly off the main JCMT beam) during the 
observations of G10.47. The CSO T225 value is often off- 
set from that recorded by the water vapour monitor, and 
does not reflect the short timescale of the opacity variabil- 
ity. The water vapour monitor, while operating during our 
observations, was new to the telescope and thus was not in- 
corporated into the data reduction pipeline. In future, the 
use of the water vapour monitor will provide more accurate 
atmospheric pwv measurements. 



13 14 15 

hours (UT) 

Figure 3. The opacity at 225 GHz measured from the CSO 
and determined through measurements of the precipitable wa- 
ter vapour by the JCMT's water vapour monitor (WVM) during 
observations of G10.47. The solid line is the CSO-measured opac- 
ity, corrected to match the JCMT observational airmass, which 
operates at fixed azimuth and updates approximately every 10-15 
minutes. The WVM is pointed slightly off the main JCMT beam, 
and updates every 6 seconds. 
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Mars observations were then used to calibrate the final, 
atmospheric transmission corrected, coadded spectra. 



3 ANALYSIS 

3.1 Source signal analysis 

We first discuss sky subtraction and instrumental contribu- 
tions to the FTS signal. We then detail calibration issues, 
which include atmospheric variation and spectral line con- 
tamination of the continuum emission. 

The measured signal S as a function of frequency v from 
an astronomical source through the FTS can be expressed 
as 

S(v) = R{v)G[na J s {v)e~ T{v) + Va J sky {v){l - e- T(v) ) 

+ (l-r)a)Jamb(v)-j'(v)] (2) 

where the contribution from any additional background has 
been neglected. Here, r(y) is the atmospheric optical depth 
as a function of frequency u, the telescope aperture efficiency 
is given by r) a , R(v) is the responsivity or response function 
of the FTS, and G is the detector gain when observing an 
astronomical source. For all observations, the gain was set to 
G = 10 4 . J a (v) is the true spectrum of the source, J s k y {v) is 
the true spectrum of the atmosphere, J a mb{v) is the power 
received from ambient temperature surfaces of the telescope, 
and J '(u) is the power received through the 2 nd port of the 
FTS, which is differenced in the interferometric measure- 
ment iNavlor et al.ll200(il) . The first term therefore repre- 
sents the power received from the astronomical source itself 
through the atmosphere and instrument, the second term 
represents the power received from the atmosphere and the 
third term represents the power received from the surfaces 
of the telescope, which are at ambient temperature. 

Observations at the background position produce solely 
the intensity received from the sky, S s k y (v), and ambient 
surfaces: 

S sky (v) = GR(v)[ria J aky {v)(l-e- T(v) ) 

+ {l-V«)Jamb(v)-j'(v)] (3) 

If the atmospheric temperature and opacity are assumed 
to remain unchanged during each pair of source and back- 
ground scans, the difference of the on- and off-source scans 
gives the observed power of the source itself, Ss(v): 

S.(v) = S{u)~S sky {u) = R{u)Gr } aJs{v)e- Ts{v) (4) 

Sky brightnesses at 345 GHz for the typical pwv levels dur- 
ing our observations of ~ 1.1mm are ~ 800 Jy /beam, an 
order of magnitude larger than our brightest source, making 
accurate background subtraction essential. During the ob- 
servations, measurements of the precipitable water vapour, 
from the CSO radiometer, and observational airmass were 
recorded. With this information, the atmospheric model de- 
scribed in Section 2 was used to determine the transmission 
of the atmosphere in the direction of the source, and thus the 
contribution of the term e -1 "'"' to the observed source in- 
tensity S 3 {y). Typical values of T225 during the observations 
were ~ 0.07, leading to atmospheric transmission levels at 
345 GHz of ~ 75%. 

Observations of Mars and a brightness temperature 
model were used for calibration of the spectra obtained by 



Table 2. Instrumental and observed S/N values 



Object 


Time on source 
(hours) 


S/N /ns 


S/N 06s 


G10.47 


3.5 


147 


8.5 


G31.41 


1.0 


99 


8.0 


G12.21 


4.0 


73 


7.4 



The instrumental and observed signal-to-noise ratio for each 
source. 

the instrument. At the time of our observations, Mars had 
an angular diameter of 8.97". The Mars observations were 
done in the same on- and off-source pattern as the HMC 
observations, and the intensity Sm,(y) received from Mars 
after sky subtraction is thus 

S m (u) = S{v) - S sky {v) = R{u)Gr/ a J m (^)e" TmM (5) 

The brightness temperature of Mars, Tg, m , on the date 
of our observations was calculated using the STARLINK 
FLUXES program JPrivett. Jenness &: Matthew J Il998h . 
The Mars flux, J m {v), at the telescope can then be cal- 
culated: 

Jm{v) = -5= (6) 

The true spectrum of the hot core sources, J a (v), can 
then be determined: 

^ = (7) 

The final 350 GHz spectrum of each source can be seen 
in Figure^Jin units of Jy /beam. The derived intensity near 
the edges of the band has larger error due to the strong 
water vapour transitions seen in Figure which define the 
350 GHz spectral window. Variations in the column abun- 
dance of water vapour between source and background mea- 
surements result in greater noise at the steepest part of the 
water vapour line profile. 

3.2 Signal-to-noise and atmospheric effects 

The instrumental signal-to-noise (S /N) ratio was calcula ted 
from the MZ FTS characteristics in lNavlor et alJ J20031) for 
the spectral resolution and specific observing times for each 
source, and is solely due to the detector noise of the FTS. 
The S/N ratios of the final spectra are compared with the 
instrumental S/N ratios in Table |5J with the average S/N 
achieved across the band a factor of order ~ 10 lower than 
the instrumental value for all sources. With all detector ef- 
fects accounted for in the instrumental S/N ratio, the sensi- 
tivity degradation is solely due to the atmospheric variance 
during our observations. 

At submillimetre wavelengths, the brightness and vari- 
ability of the sky is an overwhelming obstacle to accurate 
flux calibration, and observations generally require some 
form of chopping between source and background measure- 
ments in order to remove the sky emission. Ideally, the time 
between source and background measurements is as short 
as possible to limit the level of sky variance between ob- 
servations j Archibald et alJl2002l) . While our observational 
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Figure 4. Final calibrated spectra of the hot core sources 
(top to bottom: G10.47, G31.41, G12.21) with la error bars at 
1.5 GHz resolution. The dash-dotted line is the best fit linear 
slope. The dashed lines indicate the fitting window. Arrows in- 
dicate the positions of five molecular emission lines, listed in the 
text, expected to be strong in HMC. With the exception of the 
CO(3-2) line at 345.6 GHz, no spectral features are coincident 
with the expected bright line positions. 



method was designed to minimize the effects of the ever- 
changing atmosphere on the data by switching between 
source and background observations every 90s, the atmo- 
sphere still contributed substantially to our uncertainties 
for several reasons. First, the opacity varies on timescales 
shorter than 90s, as shown in Figure Second, delays in 
slewing the telescope to the background position resulted in 
slight differences between the airmass of the background and 
source scans, introducing a small systematic effect by effec- 
tively offsetting the opacity in the background scan relative 
to the on-source scan. 

As a quick test, we can think of a small increase in 
opacity as a small increase in the precipitable water vapour 
(with the caveat that there are other minor contributions to 
the opacity). Using atmospheric emissivities calculated for 
various pwv levels with ULTRAM, we have calculated that 
a modest jump in pwv of 0.05 mm, very reasonable when 
combining small deviations between source and background 
airmass with actual changes in the atmospheric pwv lev- 
els, increases the atmospheric emission by ~ 35 Jy/beam at 
350 GHz. The pwv values calculated from the CSO T225 for 
the observations contain jumps of this size and a few even 
larger. For this reason, stable weather is more crucial to this 
study than particularly good weather. A change in atmo- 
spheric emission of 35 Jy is comparable to the flux of our 
brightest source, G10.47, and brighter than the others (see 
Table 0. As a consequence, some of the subtractions of the 
source and background spectra produced unphysical, nega- 
tive flux measurements, while others resulted in impossibly 
high source brightnesses. When the spectra are coadded, 
however, on the whole these overly positive and negative 
scans will average out unless there is a large systematic off- 
set between source and background observations or the pwv 
levels in the atmosphere are systematically increasing or de- 
creasing. 

While there was a small systematic offset in the airmass 
between the on- and off-source scans, it was not large enough 
to produce significant changes in the sky emission except in 
the case of G31.41, which was observed while setting. For 
this source, the offset may have decreased the observed flux 
by at most 1-2 Jy/beam. During observations of G10.47, the 
recorded atmospheric opacity steadily increased, although it 
remained low. After atmospheric transmission corrections, 
the recorded intensity of this source near the start of the 
observations was greater by ~ 10% on average than that 
recorded later. Again, this may have decreased the observed 
flux by a few Jy/beam at most, and is the largest contributor 
to the uncertainty per data point in the spectrum. There do 
not appear to be any systematic trends in the atmospheric 
opacity for the G12.21 observations. 

To estimate the effects of the varying atmosphere on 
our spectra, and to determine the level of any systematic 
effects which may have been introduced into the data, we 
divided the sky-subtracted spectra into two groups for each 
source corresponding to observations made in the first and 
second halves of each night. We then coadded the spectra in 
these groups, and ratioed the coadded spectrum of the first 
half with that of the second. Figure shows this ratio for 
Mars, our calibrator, and G10.47. It is clear that while there 
is a small shift in overall flux between the early and later 
observations, within the band there is little overall change in 
the continuum slope. The fainter sources showed similar re- 
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Table 3. Comparison of source brightnesses from SCUBA and 
the FTS 



330 335 340 345 350 355 360 365 

v (GHz) 




0.6- G10.47 



Figure 5. Sky-subtracted spectra, divided into two groups cor- 
responding to observations made in the first and second halves of 
each night, then coadded and ratioed for Mars and G10.47. The 
dash-dotted line indicates the best-fit slope of the ratioed data 
for G10.47, showing that while the variable atmosphere may in- 
troduce offsets into our final observed intensity of the sources, it 
has not significantly affected the slope of the continuum emission. 



suits, and we thus expect that while the variable atmosphere 
may introduce offsets into our final observed intensity of the 
sources, it has not significantly affected the slope of the con- 
tinuum emission. 



3.3 Signal strength 

The final brightness values recorded by the FTS for each 
source are listed in Table [3] with calculated uncertainties. 
The uncertainty in the final flux value comes from the uncer- 
tainty in the observed flux of the source, S 3 (v), the observed 
flux of Mars, Smiv), the uncertainty in the Mars brightness 
temperature model, and the uncertainty in the atmospheric 
transmission, e~ T ^\ From the S/N of the observations, the 
uncertainty in S s (v) is ~ 11 — 12% for each source. We 
take an uncertainty in the atmospheric transmission of 10%, 
while the estimated uncertainty for the Mars model is of or- 
der 5%. The uncertainty in S m (i') is negligible compared to 
the other factors due to Mars' brightness. Added in quadra- 
ture, this gives an uncertainty in the calibrated source flux of 
16%. Additionally, the error lobe of the 850 urn (350 GHz) 



Object 


I v (SCUBA) 


Iv (FTS) 




Jy/15" beam 


Jy/15" beam 


G10.47 


39.0 ±5.9 


44.0 ± 8.8 


G31.41 


24.2 ±3.6 


24.6 ±4.9 


G12.21 


13.0 ±2.0 


27.2 ±5.4 



A comparison of the 850 um source brightnesses determined us- 
ing the FTS with the peak brightnesses recorded by the Submil- 
limctrc Common User Bolometer Array. The uncertainties in the 
SCUBA fluxes are 15%, while the uncertainties in the FTS fluxes 
are 21%. 

JCMT beam may contain up to ~ 5% of the source flux, 
making the final uncertainties of order 21%. 

Through calibration with Mars, the observed fluxes of 
G10.47 and G31.41 agree very well with those recorded by 
SCUBA; however, G12.21 is much brighter in our observa- 
tions than the submillimetre continuum measurements indi- 
cate, and the values are not within uncertainties. The peak 
fluxes from the archived SCUBA data for G12.21, which 
are from two separate observations, agree within uncertain- 
ties, suggesting that a variation of the source brightness is 
unlikely. G12.21 is the faintest source by the SCUBA mea- 
surements, and is thus more susceptible to atmospheric sub- 
traction problems than G10.47 and G31.41. Since the largest 
effect of the variable atmosphere on our observations is to 
introduce an offset into the overall continuum, and because 
this offset is cumulative, it is possible that the final flux for 
this source has been inflated by the bright, variable atmo- 
sphere. 



3.4 Molecular line contamination 

The sources observed in this study are hot molecular cores. 
The moderate temperatures of these objects (T > 100 K) 
cause molecules and elements frozen onto the surface of dust 
grains during the colder collapse phase to evaporate, releas- 
ing these species into the gas phase. The resulting emission 
lines are often strong in the 350 GHz spectral window for 
HMC, and can therefore contaminate measurements of the 
submillimetre continuum. The source s in this study have 
been detected in many molecular lines ijHatchell et al]ll998Fl 
from multiple species present in the 350 GHz atmospheric 
window. 

At our low spectral resolution only very bright, rela- 
tively broad lines would contribute noticeably to the dust 
continuum spectrum. Our spectra show little evidence of 
bright line contamination, as shown in Figure |1] where 
we have indicated the expected positions of 13 CO (3-2), 
CS (7-6), CO (3-2), HCN (4-3), and HCO+ (4-3) emis- 
sion lines. These are the strongest, broadest lines ob- 
served in the 350 GHz window of the hot core G5. 89-0. 39 
jThompson fc Macdonardll999Fl . which we expect to be sim- 
ilar to the cores studied here. There is some evidence of 
small peaks in all the spectra, one of which appears to be at 
the CO(3-2) transition frequency (345.6GHz); however the 
data do not allow a conclusive identification. None of the 
other observed peaks in the spectra are coincident with ex- 
pected bright line positions. Overall the peaks are small and 
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do not significantly influence the continuum slope, and the 
contamination of the continuum spectrum by single bright 
emission lines is therefore not a large problem for this anal- 
ysis. The concern for this study is then the contamination of 
the continuum spectrum by a significant number of low-level 
emission lines that may be present and are not immediately 
visible, but which may influence the slope of the continuum. 

The frequency coverage of molecular line studies for 
these sources is much less than that covered by SCUBA and 
the FTS, so no absolut e estimates of the lev el of contami- 
nation have been made jHatchell et aljfeoOCfl . Observing at 
350 GHz with the JCMT, the conversion factor from bright- 
ness temperature T, in Kelv in to flux S in Jansky s is approx- 
imately S(Jy) ~ 13T b (K) faatchell et alJll998T) . A bright- 
ness temperature of 40 K thus corresponds to ~ 500 Jy. The 
line width of strong molecular features, however, is seldom 
greater than 15kms -1 (at least for the bright core of the 
line) while the spectral resolution of the FTS for this study is 
~ 1500 km s~ , so that individual lines are diluted by at least 
a factor of 100. Strong molecular lines, with peak brightness 
~ 40 K, will thus add approximately 5Jy to the measured 
flux in a given spectral bin. Therefore, we do not expect con- 
tamination of the spectra by individual lines by more than 
approximately 10-20%. 

The contamination of the submillimetre continuum 
of hot core sources by molecular lines can range be- 
tween 10% and (in extreme cases) 60% of t he to - 
tal integrated flux (iGroesbeck. Phillips fc Blakd Il994l) . 
Ijohnstone. Boonman &^an^ishoedTffiOo3h ~t^d that the 
350 GHz continuum emission from their observed protostel- 
lar sources, however, is never dominated by molecules other 
than CO (although line emission from HCN, HNC, CN and 
CH3OH can contribute a substantial fraction of the con- 
tamination for more energetic sources) and line contamina- 
tion in general is typically less than ~ 10% at 350 GHz, 
even for photon-dominated regions. The hot core regions of 
the sources observed in this study amount to a small frac- 
tion of the 350 GHz beam (1-2" in diameter compared with 
15" FWHM), whereas the more extended emission is from 
cooler material in which fewer molecules will be excited. 
Overall, we argue that molecular line contamination of our 
low resolution 350 GHz continuum is minimal, and does not 
significantly affect the analysis. We recognize, however, that 
it is important to consider the influence of these low-level 
molecular lines when making any determination of the slope 
of the dust continuum emission using low spectral resolution 
submillimetre observations. 



4 SPECTRAL INDEX RESULTS 

4.1 Calculation of the spectral index 

From Equation Q the dust emission S(v) oc v 1 , where 
7 = p + a is the sum of the dust spectral index /3 and 
the frequency dependence a from the Planck function. For 
small ranges in frequency A v around a central frequency vq, 
we can expand this relationship: 

S(v Q + Av) <x (md + A u) 1 

ex ^(l + j^) (8) 



Table 4. Spectral indices 



Source 


"YFTS 





1SCUBA 


G10.47 


3.6 ± 1.2 


1.6 ± 1.2 


3.4 ±0.6 


G31.41 


4.1 ± 1.2 


2.1 ± 1.2 


3.3 ± 0.6 


G12.21 


3.3 ± 1.1 


1.3 ± 1.1 


4.1 ± 0.6 



Thus, for small Av /uo, the dust emission is expected to 
increase linearly with frequency. The FTS passband at 
350 GHz is ~ 30 GHz wide, or ~ 10% of the observing fre- 
quency, making the approximation valid. In practice, this 
window was made slightly smaller due to a decrease in S /N 
caused by increased levels of contaminating atmospheric flux 
near the edges of the band, leaving A v/u ~ 0.08 on average. 
A x 2 minimizing linear fit routine was used to determine the 
slope of the continuum emission inside the passband. All re- 
sults quoted here have been determined from the unapodized 
data to retain the highest spectral resolution. The results 
varied slightly depending on the fraction of the band used 
(due to increasing error bars at the band edges when using 
more of the band, and conversely due to a smaller lever arm 
when using less of the band); however any variations were 
within uncertainties. 

Column 2 of Table^Jshows the best fit continuum emis- 
sion slope 7 for each source with uncertainties determined 
through the \ 2 fitting routine, and Figure [I] shows the final 
spectrum for each source overlaid with the best linear fit. 
The average reduced x 2 value for the fits was ~ 0.6, indicat- 
ing a good fit to the data within the band. We expect, for a 
good linear fit, that the reduced \ 2 should be close to unity 
for a sufficiently large number of data points. The small val- 
ues of the reduced \ 2 m our fits suggest that the uncertainty 
per data point has been overestimated by ~ 20% for each 
source. Since the uncertainty in the calculated slope varies 
as the uncertainty per data point, this implies that the un- 
certainties in our slopes may be overestimated by ~ 20%. 

For emission in the Rayleigh- Jeans limit, a = 2 and 
thus the dust spectral index f3 = 7 — 2. In order for the 
R-J approximation to apply at 350 GHz, hv/kTd <C 1, so 
T d > hv/k ~ 17 K. Although HMC are warm, they are 
surrounded by a shell of cold dust. The observed emission 
is integrated over the line of sight and the beam, which at 
15" FWHM is much larger than the hot cores. There is there- 
fore substantial contribution to the observed emission from 
the colder outer regions along the line of sight. Models of 
the emission from these sources which take the source tem- 
perature and density profiles into account, however, appear 
to produce spectra which are in the R-J limit at 3 50 GHz 
jHatchell et al.ll2000l : iHatchell fc van der TaJdl2003ft . When 
considering only the peak flux, a greater proportion of the 
emission from the warm inner regions is contained within 
the single beam, and the R-J approximation still applies. 
Column 3 of Table 2] lists the dust spectral index /3 calcu- 
lated for each source assuming the 350 GHz emission is in 
the R-J limit. We find in the range of 1.3 — 2.1 with a 
mean of 1.6. 
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4.2 Comparison with SCUBA data 

As discussed in Section 1, the spectral index of dust emis- 
sion can also be calculated using observations at two widely 
separated wavelengths. From Equation^ we see that we can 
solve for 7 = /3 + a with a ratio of the flux at two different 
wavelengths: 



7 = 



Hga / Si 
\n(i/2 I v\) 



= /3 + a 



(9) 



For our sources, archived SCUBA data at 450 /J.m (660 GHz) 
and 850 /im (350 GHz) are publicly available. In order to 
properly compare the flux at the two wavelengths, the 
8.5" resolution 450 [im data were first convolved to the 
15" resolution of the 850 u.m data. The ratio of the peak 
brightnesses was then used to calculate 7 from Equation 
Column 4 of Table 0] gives the 7 value, 7 scuba, calcu- 
lated from the ratio of the 450 u.m and 850 u.m fluxes for 
each source. Within uncertainties, "/scuba is consistent with 
7 fts • The uncertainties in this calculation are dominated by 
the uncertainty in the source flux at both wavelengths. 

If the R-J approximation applies at both frequencies 
v\ and U2, then a — 2 and /3 = "/scuba — 2. As shown 
previously, the emission at 350 GHz is likely in the R-J 
limit, but for the R-J approximation to apply at 660 GHz, 
Td S> hv/k ~ 32 K. For the average dust temperatures 
likely associated with these hot molecular core sources, the 
660 GHz emission cannot be assumed to be in the R-J limit. 
In this case, less power is contributed to the dust emission 
from the Planck function, and a < 2. Using this informa- 
tion, we can use the 350 GHz (850 u.m) slope determined 
using the FTS and "/scuba to estimate the temperatures of 
the observed HMC. 

The spectral indices calculated for G10.47 agree well 
with a ~ 2, suggesting that the average dust temperature 
for this source is high enough (Td > 60 K) for the emission at 
both wavelengths to be in the R-J limit. G12.21 is also likely 
warm, while the results for G31.41 indicate a lower average 
dust temperature. The uncertainties in these temperature 
estimates are very large due to the large uncertainties in 
the evaluated slopes; it is clear, however, that with smaller 
uncertainties in the value of /3 determined from FTS obser- 
vations, it will be possible to independently determine both 
Td and /3 of the dust opacity in star forming regions. 



4.3 Interpretation 

Our results are consistent with the values of f3 
found in embedded UCHII regions and other HMC. 
Fitting models of H MC to observational data have 
found = 1.6-2 JOsorio. Lizano fc D'Alessid Il999t 
IChurchwell. Wolfire fc Wood! Il99fll . In a study of 17 
UC H II regions, iHunten (I1998T) finds an average spectral 
index of f3 = 2.0 ± 0.25 through greybody modelling of 
IRAS, submillimetre and millimetre observations. Observa- 
tions and modell ing of six compact HII regio n s, one of which 
was G10.47, by iHoare. Roche fc Glencrossl ( 1991) find an 
average /3 ~ 1.5. They note, however, that the frequency de- 
pendence is somewhat affected by the assumed density dis- 
tribution of the model. The dense ridge in Orion has 13 ~ 1.9 
near embedded infrared sources llGoldsmith. Bergin fc Lia 
Il997l) . Overall, the spectral indices of these warmer, likely 



more evolved star forming regions tend to be higher than 
those found in cooler, less evolved dense cores, indicating 
that dust properties have changed through some mechanism 
which appears to correlate with protostellar age. 

Several dust models predict a dust emissivity frequency 
dependence similar to that found here, but our observations 
are not sensitive enough to determine a most likely dust 
composition. Higher values of /3 have been associated with 
dust grains covered in thick ice mantles iAanestadl Il975l : 
iLis fc Mentenl [l998'l . while many dust grain models with 
combinations of silicate and graphite compositions find a 
spectral index of /3 ~ 1.5 — 2. lOssenkopf fc Hennind <ll994f) 
tabulate the dust opacity properties for nine different dust 
distributions, including bare grains versus grains with ice 
mantles (both thick and thin), and various degrees of dust 
coagulation. The emissivity results presented in this paper 
are unable to definitively rule out any of the models. The 
best fit, however, is for bare grains in relatively dense envi- 
ronments, w here coagulation is mor e pron ounced (column 2 
of Table 1 in lOssenkopf fc Henningl lll994h ). 

While we had hoped to determine 7 to an accuracy 
of ±0.1, the variable atmosphere at 350 GHz degraded the 
achieved signal-to-noise ratio by a factor > 10 compared 
with the instrumental S/N, leading to uncertainties in our 
calculation of j3 of approximately ±1.2. Even with the at- 
mospheric complications, however, the uncertainties deter- 
mined here are of the same order as those found when calcu- 
lating the dust spectral index through multiple wavelength 
observations of single sources. These results show that with 
better sky subtraction techniques, it should be possible to 
determine the spectral index of bright astronomical sources 
with high accuracy using an FTS. Future instruments with 
greater sensitivity and precise sky subtraction will enable 
the measurement of spectral indices of much fainter sources. 
These results show that careful sky subtraction is essential 
to future dust continuum emission studies with an FTS. 



4.4 An imaging FTS for SCUBA-2 

SCUBA-2 is the next generation submillimetre instru- 
ment currently in de velopment to replace SCUBA in 2006 
^Holland et alJ 12003ft . An imaging FTS (IFTS) is be- 
ing developed as an ancillary instrument to SCUBA-2 
jNavlor fc Goml I2004T) . Based on detector noise calcula- 
tions alone, the IFTS will be at least 10 times more sen- 
sitive per pixel th an the single bolome ter FTS used in 
these observations jGom fc NavloJ 12004ft . The extremely 
accurate sky subtraction possible with the IFTS will pro- 
vide an even more substantial increase in the observational 
precision. When also considering the variable atmosphere 
noise, the per-pixel sensitivity of the IFTS increases to 
at least 100 times greater than in the observations pre- 
sented in this paper . Using the detector sensitivities in 
iGom fc NavloJ (2004), we have calculated that it will be 
possible to determine in 12 hours of observing with the 
IFTS the spectral index 7 = /3 + a of the submillime- 
tre dust emission at 850 [im to ±0.1 for sources of only 
400mJy/beam in brightness, and to ±0.3 for sources only 
~ 140mJy/beam in brightness. Similar calculations for ob- 
servations at 450 n.m are listed with those at 850 urn in Table 
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Table 5. Source brightnesses required to determine 7 in 12 hours of observing 



Wavelength 


Frequency 


Resolution 


12 hrs la AT 


Source brightness 


Source brightness 


yim 


GHz 


GHz 


mK 


7 ±0.1 (Jy/beam) 


7 ± 0.3 (Jy/beam) 


450 


660 


3.0 


1.44 


4.7 


1.6 


850 


350 


3.0 


0.2 


0.40 


0.14 



Source brightnesses required to determine the dust spectral index 7 in 12 hours of observing using an 
IFTS on SCUBA-2. 



5 SUMMARY 

We have determined the spectral index of the dust emis- 
sion of three hot molecular cores solely within the 350 GHz 
(850 |im) passband using an FTS on the JCMT. We find an 
average dust spectral index (3 ~ 1.6, in agreement with the 
spectral indices calculated using archived SCUBA data at 
450 U-m and 850 um. These results for /3 are consistent with 
measurements of the spectral index in other HMC. The un- 
certainties in p from these observations are dominated by 
the difficulties in subtracting the very bright and variable 
submillimetre atmosphere from the data. With better sky 
subtraction techniques, as will be possible with an imaging 
FTS on SCUBA-2, these uncertainties will be greatly re- 
duced. The per-pixel sensitivity of the imaging FTS planned 
for SCUBA-2 will be greater by a factor of ~ 100 than that 
of the single pixel detector used in this study, and will allow 
the determination of the spectral index of dust emission of 
sources with brightnesses of only a few hundred mJy/JCMT 
beam. 
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